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NECESSARY CONDITIONS FOR 
OPTIMAL LUNAR TRAJECTORIES WITH 
DISCONTINUOUS STATE VARIABLES AND 
INTERMEDIATE POINT CONSTRAINTS 

by 

W. E. Miner* and J. F. Andru$ + 


SUMMARY 

The guidance regime for an optimal multi-stage lunar 
trajectory is derived by applying the mathematics of the cal- 
culus of variations as established by Denbow for the genera- 
lized problem of Bolza. The steering angle programs for 
four constant thrust level phases and the times to initiate 
and terminate the two coast phases are determined in order 
to place maximum payload into a specified lunar orbit. The 
problem considered here is to determine an optimal trajec- 
tory consisting of six sub-arcs utilizing three vehicle stages 
on which maximum payload is transported from an exo-at- 
mospheric point near the Earth to a specified lunar orbit. 

The intermediate point constraints include two points at which 
stages are separated and mass discontinuities occur, an 
Earth parking orbit of specified energy and angular momen- 
tum magnitude, and four thrust magnitude levels. The 
Euler- Lagrange equations determine the optimal steering for 
the thrusting phases and the Denbow transversality equations 
are used to calculate the discontinuities at the ends of the sub- 
arcs. This method is applied here and the equations neces- 
sary to solve this problem using a high-speed computer are 
derived. 


* Chief, Guidance Theory and Trajectory Analysis Branch, Electronics Research 
Center, NASA, Cambridge, Mass. Formerly with the Aero-Astrodynamics 
Laboratory, Marshall Space Flight Center, NASA,. where the work represented in 
this paper was initiated. 

+ Northrop Space Laboratories, Huntsville, Ala. Work performed while with the 
Huntsville Operation of the General Electric Computer Department under contract 
to the Computation Laboratory, Marshall Space Flight Center, NASA. 


INTRODUCTION 


The calculus of variations is a well established mathematical tool for the deter- 
mination of control functions and control parameters for optimal trajectories. The 
Euler- Lagrange equations yield the necessary conditions for extremal solutions to the 
split boundary value problem of Bolza. Extensions to the theory to provide for inter- 
mediate point constraints, such as those introduced into the system by the separation 
of vehicle stages or the initiation of a coast phase, were made by Denbow. 

In general, the path of an aerospace vehicle can be described by a set of n non- 
linear differential equations 

x± = f i (x(t), a (t ) , a) i = 1,2,..., n 

where 

x (t) is an n vector of state variables, 

(x L , x 2 , . . . , x n ) 
a (t) is an m vector of control functions, 

( a l ' a 2 ' *•*' °‘m ) 

IT is a k vector of control parameters, 

i \ T 

(a^ , a 2 , ■ • • i 

and T means transposed. A particular path results when the initial conditions x(t Q ), 

the control functions a(t) for t < t < t c , and the control numbers, a, are given 

o— r 

where t Q and tj are the initial and final times. 

An optimal trajectory is one which results when a (t) and a are chosen so that 
a specified terminal quantity $ (x (t , t^) , called the pay-off, is maximized or 
minimized and a specified p vector of terminal constraints 

(x (t f ) , t f ) =*0 j = 1,2, ...,p 

is satisfied. 
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I 


The calculus of variations provides for the solution of such problems by intro- 
ducing the n vector function X (t) of Lagrange multipliers, which satisfies the set 
of linear differential equations 


n 



subject to the m relations 


n 



as well as appropriate terminal conditions. 

The problem is complicated by the presence of intermediate point constraints. 
These divide the trajectory into a sequence of sub-arcs. Specifying the values of 
state variables or control functions at intermediate points imposes constraints on the 
Lagrange multipliers at these points. The calculation of these constraints at the ends 
of the sub-arcs so that the composite trajectory is an extremal was made by Denbow. 
His results can be applied to intermediate points at which state variables or control 
functions are subject to specified finite discontinuities, e.g. , the cases of vehicle 
staging and switching from powered flight to coast. These necessary conditions at the 
ends of the sub-arcs, which are called transversality conditions, can be written in an 
elegant form as upper bounds on the rank of a prescribed matrix. 

This paper is an application of the theory of optimal control with intermediate 
point constraints to a specific lunar trajectory problem. This multistage trajectory 
is reduced to a two-point boundary value problem and the results comprise the 
necessary conditions for an extremum. In order to obtain a specific optimal lunar 
trajectory, it is now necessary to solve the two-point boundary value problem. 
Numerical methods for converging to a solution, i.e. , an extremal trajectory, have 
been described widely in the literature. 


1 Denbow, C. H. : 


A Generalized Form of the Problem of Bolza. Contributions to the 
Calculus of Variations, University of Chicago Press, Oct. 1937. 
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A list of symbols used in this Technical Memorandum is presented below: 


A 


A* 1 
A** r 



F i 


G 

H i 

H* 

I 

’sp 

M 

M* 

O 

R 

S 


List of Symbols 

Augmented matrix shown in Figure 3 

Sub-matrices of A defined on pp. 24, 25, and 27 

Matrix which relates Earth-centered and moon-centered coordinate 
systems 

Sub-matrix of A associated with the parking-orbit constraints 
defined on p. 23 

Non- singular submatrix of C 3 

Sub-matrix of A associated with the lunar orbit constraints 
Non- singular sub-matrix of C 7 

Matrices defined on pp. 32 and 33 

Matrix defined on p. 35 

Vehicle thrust magnitude during i tfl flight phase 
i = 1 , 2 , ...,6 

Gravitational force 

The Hamiltonian evaluated at point i 

Expression defined on p. 33 

Identity matrix 

Specific impulse 

Sub-matrix of A defined on p. 28 
Product matrix defined on p. 33 
Matrix of zeros, dimensions indicated 
Distance to entry point of parking orbit 
Euler- Lagrange sum defined on p. 15 



Flight duration from exit of Earth parking orbit to entrance of lunar 
orbit or matrix transpose when used as a superscript 

Speed at entry point of parking orbit 

Cartesian coordinate components of G 

k vector of control parameters 
General element of A 

Quantities constrained at sub-arc boundaries 
Eccentricity of lunar orbit 
n vector of rates of state variables 
Gravitational constant at surface of Earth 
Semi-latus rectum of lunar orbit 
Variable vehicle mass at point i 

Magnitude of position vector in moon-centered coordinate system 
Time and state variables at sub-arc boundaries 
Time at point i 

Vehicle velocity components in Cartesian coordinate system 
Vehicle velocity vector at point i 

Vehicle velocity vector in moon-centered coordinate system 

Vehicle position components in Cartesian coordinate system 

General n vector of state variables 
Vehicle position vector at point i 

Vehicle position vector in moon-centered coordinate system 

Position of moon with respect to Earth-centered reference system 
at point 7 
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$ 

X, 


♦i 


Quantity defined on p. 17 
m vector of control functions 
Fuel burning rate during the i*"* 1 flight phase 
+1 , used to indicate sign ambiguity 

The i m Lagrange multiplier 

T 

The vector ( X ^ , X ^ , X ^ ) evaluated at point i 

Gravitational constant of moon 

T 

The vector ( X ^ , X , X g ) evaluated at point i 
Expression defined on p. 33 

Expressions defined on p. 36 

General variable used to represent, in turn, each of the state 
variables and control functions 

Terminal quantity to be maximized or minimized, the payoff 

Thrust angle in pitch\ control variables 

Thrust angle in yaw I shown in Figure 2 

General terminal constraint 

Constraint describing specified lunar orbit 

T 

Vector • • • ^5] 



THE PHYSICAL PROBLEM 


The vehicle is a variable mass particle which moves from a point outside the 
Earth's atmosphere into a specified lunar orbit under the action of its own thrust and 
a gravitational field computed using an ephemeris. The problem is to program the 
thrust steering for the powered phases and the times of initiation and termination of 
the coast phases so that the payload placed in orbit about the moon is maximum. 

A non-rotating, Earth-centered Cartesian coordinate system is the principal 
frame of reference. The terminal orbit is given in a moon-centered system which is 
described when it is introduced. 

Conditions are imposed at seven points along the trajectory as shown in Figure 1. 
These divide the flight path into six sub-arcs on each of which the magnitude of the 
thrust is constant and designated F^, F 2 , F^, and F g on the powered arcs. The two 
coast periods correspond to F 3 = F 5 = 0. The mass rate is a constant on each sub- 
arc. The vehicle has three stages so that the mass has a specified discontinuity at 
two points equal to the mass of the separated stage. Further, the third sub-arc is a 
parking orbit about the Earth which is only partially specified so that the remaining 
quantities may also be optimized. 



Figure 1. - Lunar trajectory for three- stage vehicle with powered and coasting arcs 
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The equations of motion are: 


u = - | sin Xp cos X y + X 

• F 

V = — COS X COS X + Y 
m p A y 

• F 

w = - sm x + Z 
m y 

x = u 
y = v 
z = w 
m = - 6 

— T _ T 

where x = (x, y, z) is the position vector, u = (u, v, w) is the velocity vector, 

F and 6 are specified constants different on each sub-arc, m is the mass, X p 

and x y are the steering angles for the thrust vector, as shown in Figure 2, andx, Y, 

Z are the components of the gravitational field. We seek x (t) and X y (t) and the 

times to initiate and terminate the coast phases. ** 


z 



Figure 2. - Definition of v and v 

P *y 
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The intermediate point and boundary constraints are specified at seven points, 
designated point 1, point 2, . . . , point 7, which divide the trajectory into six sub-arcs 
called First Stage, Second Stage A, Parking Orbit Coast, Second Stage B, Trans- 
Lunar Coast, and Third Stage. All quantities with finite discontinuities will exhibit 
this by plus or minus superscripts on their subscripts; e. g. , m 2 _ and m 2+ indicate 
the mass at point 2 just prior to, and just after, separation of the first stage, 
respectively. The description on conditions specified on the six sub-arcs are as 
follows. 


First Stage 


Initially the vehicle is located at point 1, outside the atmosphere of the Earth 
with initial time, position vector, velocity vector, mass, thrust magnitude, and 
burning rate specified. Thus we have t^, Xp u^ m^ F^, and 8^ given. This 
stage terminates at a given time t 2 when the empty fuel tank is separated. There- 
fore, the mass at the end of this arc is: 

m 2 _= (t 2 -ti) + m 1 • 

Throughout this arc, F^ and are held fixed and the thrust direction is optimized. 

Second Stage A 


At point 2, discontinuities in mass, thrust level, and burning rate occur. Thus 
we have: 


m. = + Airu 

2 + 2 - 2 


2+ 


F 1 + AF 2 


2 + 


Bl + Ag 2 


where Am 2 , AF 2 , and A8 2 are given. On this sub-arc the vehicle continues to 
thrust, varying its thrust direction angles optimally until it arrives in a parking orbit 
about the earth at point 3. The parking orbit is partially specified by requiring: 
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- 2 
x 3 . x 3 = R 



u 3 • x 3 ■ 0 

where R and V are given quantities. This is equivalent to specifying the energy 
and magnitude of the angular momentum of the orbit and that the vehicle enter the 
orbit at an apsis. The plane of the orbit and the inclination of the orbit in the plane 
are left free to be optimized. 


Parking Orbit Coast 


When the orbit is achieved at point 3, the engines are switched off and discon- 
tinuities occur in thrust level and burning rate. Thus: 

F 3 - P 2+ -F 2+ “ 0 


B 3 B 2+ -B 2+ 0 


The vehicle coasts in this orbit until such time as it is best to resume thrusting again 
in order to enter a trans-lunar coasting arc. The mass is constant until point 4 is 
reached and coasting terminates. 


Second Stage B 


Engines are turned on again at point 4 and powered flight resumes. There are 
then discontinuities in F and 6 of magnitudes F^ + and 3^ + , respectively, and 
these quantities remain constant throughout this phase. The time t 4 is a quantity to 
be optimized. There is a restriction on the difference t^ - t^ where t^ is final 
time. This difference is the time of flight from the parking orbit to the lunar orbit 
and is a given number T. Thrusting continues under optimal steering until point 5 is 
reached at which point the trans-lunar coast begins. The time t 5 is also optimized. 
At this point the mass is given by 

m 5- = m 2+ " B 2+( t 3~ t 2) - B 4+( t 5- t 4)‘ 
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Trans— Lunar Coast 


At point 5, the second stage is separated and coasting resumes. Thus: 
m 5+ = m 5- + Am 5 



e 5 = 


0 


where Anricj is specified. Amg can be considered the mass of the empty second 
stage. Any unburned fuel is treated as part of the payload. This phase takes the 
vehicle to the vicinity of the moon and terminates at point 6, optimizing the time tg. 

Third Stage 

At point 6, thrusting resumes under optimal steering. Thus, there are dis- 
continuities in F and 6 of magnitudes Fg + and 6g + , respectively, and these 
quantities are constant until the vehicle arrives in the specified lunar orbit. Point 7 
is the terminal point at which five constraints describing the given lunar orbit must 
be satisfied. Thus: 

u ~j / x ) = 0 for i = 1 , 2 , . . . , 5 . 

These equations specify a conic section expressed relative to a rotating moon-centered 
coordinate system rigidly attached to the moon. These conditions are time-dependent 
because the coordinate system rotates with the moon. Recall that, as stated in 
Second Stage B above, the time t^ was constrained to satisfy the equation 



where T is a given number. 

In summary, the purpose of this study, and the definition of optimal trajectory 
as used herein, is to determine the steering histories and switching times which yield 
a trajectory that is a mathematical extremum for placing maximum payload into a 
specified lunar orbit while satisfying the intermediate point and boundary constraints 
described above. 
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EULER-LAGRANGE EQUATIONS 


The work of Denbow has established that a necessary condition for an extremum 
of the generalized variational problem is that the classical Euler- Lagrange equations 
are satisfied on each sub-arc of a path subject to intermediate point constraints. 

Note, that to be consistent with Ref. 1, maximizing m ?) the mass of the vehicle at 
point 7, is equivalent to minimizing -m ? . 


Recalling the equations of motion in the previous section, the Euler- Lagrange 
sum is: 

S = X x (u + l sin x p cos Xy -X) 

+X 2 " I cos x p cos x y " Y) 

+X 3 (w - | sin Xy - Z) 

+X 4 (x - u) 

+X 5 (y - v) 

+Xg (z - w) 

+X 7 (m + 6) 


The Euler- Lagrange equations can be determined from the conditions stated in the 
Introduction or equivalently from: 



where x = u, v, w, x, y, z, m, X p , andX y . The resulting equations are: 
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X 1 = - X 4 


X 2 " - X 5 
X 3 = ~ X 6 


• 

x 


9X 

X 

9 Y 

- X p. 

4 

1 

9 x 

2 

9x 

3 9x 

• 


ax _ 


9 Y 

i 32 

X 5 

= " X 1 

3y 

X 2 

9 y 

“ 3 ?y 

• 


ax 


9 Y 

, az 

X 6 

= - X 1 

9 z 

X 2 

3 z 

” X 3 157 

• 

x„ 

= * I 

-X. 

sin 

X 

cos x + X cos x cos x„ + X sin x„ 

7 

2 l 

m \ 

1 


P 

y 2 p y 3 y 


- E (a, 

m \ 1 

= — / -X si 
m \ 1 


cos x cos X„ + X 0 sin x_ cos x, 


sin Xp sin x y + X 2 cos x p sin x y -* 3 cos x y 


The last two equations yield: 


tan Xp ■ - q 

X 2^3 

cos x tan x = - t 

P Y \ 2 + x 2 

A 1 + A 2 


The mathematical possibilities cos X y = 0 and F = 0 lead to trivial cases. The 
above equations can be written: 


sin x 

P 



X 


2 

2 
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where <5p = + 1 and 6y = +1. The ambiguities in the signs will be resolved in 
the next section. Note that we have the optimal steering angles, albeit in terms of 
the Lagrange multipliers. 

If we define the quantity A as: 



and 
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then the system of differential equations for the optimal trajectory can be summarized: 


F . A . 

m S y A + G 


m = -f$ 


3G 


X? - 


8x 

F 


m 


<S A 

2 y 


where the optimal steering is given by 


-1 


= tan 


-1 


X.. = sin 


6 A- 
y3i 
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WEIERSTRASS NECESSARY CONDITION 


In Ref. 1, it is proved that the classical Weierstrass condition must hold on 
every sub-arc. Applied to this particular problem, this requirement reduces to 
specif ying that the negative of the Euler- Lagrange sum, S, is a maximum with 
respect to the control functions on the optimal path. Considering only the terms 
involving Xp and Xy and replacing these by their optimal values as found in the 
last section; the quantity to be maximized is: 



This constraint resolves the ambiguity in 6 . by i S now seen to be +1. Since X y 

has been restricted to the first and fourth quadrants, the equation for cos Xy 
in the preceding section indicates that 5 y and ^ p must have the same sign. 
Therefore, ^ p = +1. The optimal steering functions are now 



where the quadrant of x p and Xy can be determined from the signs of a^> *2’ 

and X , in the usual manner with - i < y < - 

3 2 — A v — 2 
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TRANSVERSAUTY CONDITIONS 


Since the specified point constraints are insufficient to determine all required 
boundary conditions for the sub-arcs of the extremum, it is necessary to employ 
transversality conditions. We will apply the transversality relations of Denbow’s 1 
generalized problem of Bolza to this particular problem. Each of these conditions 
applies to an end point of a sub-arc. 

It can be shown that there exists a matrix, called the augmented matrix A, ^ 
such that placing the proper restrictions on its rank yields all the necessary trans- 
versality conditions. The matrix A for this problem is a ( 20 x 56) rectangular 
array with rank less than 20. Thus all (20 x 20) submatrices are singular. The 
term, a^ , i^ 19, 20 of A is a partial derivative of a point constraint with respect 
to a state variable. This matrix is shown in Figure 3. If we designate by C| the 
quantities constrained at the end points of the sub-arcs, i. e. , at point 1, point 2, . . . 
point 7, then there are 18 such terms and as i takes on the integer values from 1 to 
18, C| takes on the values tp up vp Wp xp yp zp mp t 2 , U 3 *U 3 , x 3 -x 3 , u 3 -X 3 , 
t^-t^, 4>p 2 , 4 * 4 , ifijj. Let Sj designate an element of the time or state at each 

of the critical points, i. e. , s-^ = tp s 2 = Up s 3 = Vp s 4 = Wp s 5 = Xp 

s 6 = yp s 7 = z i» s 8 = m l* s 9 = ^2’ s 10 “ u 2’ • • * s 55 = z 7’ s 56 = m 7* 
There are 8 terms at each point and since there are 7 points, there are 56 values 

for Sj . The general term in the matrix can then be defined: 


a 


ij 



i = 1,2, ,18. 

] = 1,2,. ..,56. 


The 19th row is associated with the payoff quantity, maximum terminal mass, and 
contains -1 in the last column and zeros elsewhere. The 20th row consists of changes 
in the Hamiltonian and in the Lagrange multipliers evaluated at the seven endpoints 
of the sub-arcs. It is useful to divide the seven Lagrange multipliers into three parts 
and call them: 

t The work was done in 1963 by R. W. Hunt, Southern Illinois University, in an 

unpublished report entitled "A Generalized Bolza-Mayer Problem with Discontinuous 
Solutions and Variable Intermediate Points, " for the Astrodynamic and Guidance 
Theory Division of the Aero-Astrodynamics Laboratory, George C. Marshall 
Space Flight Center, Huntsville, Ala. The paper by Denbow cited above gives 
sufficient information to construct the required matrix. 
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Figure 3. 






{X 1' X 2' X 3 ) 


( r P 2 / Pg ) ( ^ ^ / X 5 / X g) 


A with the vector u, p with x and X 7 with mass. The Hamiltonian is 

* , . 

H = -X. • u - P* x +A^B • 

The 56 quantities in the 20th row are grouped into sets of eight, each set evaluated at 

one of the seven critical points. Thus the first eight elements are -H,— A 1 ,-p , -A_ 

~ T — ffi / 

each evaluated at point 1, the next eight elements are -AH, -(A A ) , — (Ap). and 

-AA 7 evaluated at point 2 and similarly for points 3, 4, 5, and 6. The last eight 
elements are H, A T , p T , A^ evaluated at point 7. 

We now impose upon A the condition that its rank be less than 20. By setting 
the determinants of all 20 x 20 submatrices of A to zero, we will have all the 
necessary transversality conditions. Many of these determinants will be trivially 
zero and provide us with no information. 



associating 
given by 


A careful study of the augmented matrix A reveals that many of the submatrices 
belong to a special form that makes obtaining some of the transversality conditions a 
matter of inspection. Consider the (3x6) submatrix associated with the parking 
orbit constraints, 



2a 3 

2v 3 

2W 3 

0 

0 

0 

C 3 S 

0 

0 

0 

2x 3 

to 

*< 

OJ 

2z 3 


X 3 

^3 

Z 3 

U 3 

v 3 

W 3 


This matrix must have rank 3 for these three constraints to be independent. Let Cg* 
be (3x3) non-singular sub-matrix of Cg. In the same way, let the (5x7) sub- 
matrix C rj be that one associated with the terminal lunar orbit at point 7. These 
represent five independent constraints and so must have rank 5. Let C* be any 
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(5x5) non-singular submatrix of C7. We can now write a very convenient (19x19) 
submatrix of A using all rows except the last one. Thus: 


A* = 
(19x19) 


I 

(8x8) 


~i r 

(lxl). 


0 

(8x11) 

0 

(1x10) 


0 

(11x8) 


c * 

3 

(3x3) [ 


0 

(10x1) 


-1 

(lxl) 


1 ° 1 

| ( 7x3 ) 


C * 

f 

(5x5) 


0 

(6x1) 


0 

(3x7) 


0 

(1x6) 

0 

(5x1) 

-1 

(lxl) 


0 | 
(1x5) 


Note that to obtain the (1x6) 0 matrix, the 25th column was added to the 49th 
column, an operation which does not change rank. 

Since all (20 x 20) submatrices of A are singular, we can augment A* with 
any (20 x 1) column not already partly used in A* and the corresponding part of the 
last row, obtaining a singular (20 x 20) matrix. In particular, if this (20 x 1 ) 
column contains zeros in the first 19 rows, the last element in that column must 
vanish since A* is non-singular. There are 32 such columns, and thus we can 
obtain, by inspection, 32 transversality conditions: 

( AX ) ^ = 0, ( Ap ) ^ = 0, ( AX ^ ) ^ = 0 for i = 2, 4, 5,6. 

(AH) i =0 for i = 3,5,6 

( A X ? ) 3 = 0 
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This means that all the Lagrange multipliers are continuous at points 2, 4, 5, 
and 6; H is continuous at points 3, 5, and 6; and is continuous all along the 
trajectory. 

Consider now the (16 x 16) sub-matrix: 


A** = 
(16x16) 



I 

— 

I 

1 

1 

0 

(8x8) 

1 (8x8) 


1 1 

i ° 


, (lxl) 
1 

1 (1x7) 


T -i ! o 

0 

1 

, (lxl) (1x6) 

(8x8) 

1 A 

1 ' 1 n 


1 0 

1 c 7* - 0 


j (7x1) 

| 0 1 (5x5) 1 (5x1) 


1 

| ( 6xl > | o -1 


1 

| (1x5) 1 (lxl) 

1 1 _ 


We need to form all possible (20 x 20) matrices using A**, contributions from the 
(3x6) matrix of parking orbit constraints Cg and the last row. Since A** is non- 
singular, the (4x6) matrix: 


2u 3 

2v 3 

2W 3 

0 

0 

0 

0 

0 

0 

2x 3 

2y 3 

2Z 3 

X 3 

y 3 

Z 3 

U 3 

V 3 

W 3 


(-ax 2 ) 3 

(-ax 3 ) 3 

(~ Ap l) 3 

(-Ap 2 ) 

(~Ap 3 ) 


has rank less than 4. Thus, all (4x4) determinants vanish. This constraint 
yields three transversality conditions. We can write the above (4x6) matrix in a 
more concise manner as: 
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25 3 

0 

(1x3) 

(1x3) 

0 

2x 

(1x3) 

3 


(1x3) 

T~ ~ ~ ~ 

3 

U 3 

(1x3) 

(1x3) 

<-Ix ) 3 

1 (- Ap ) 3 

(1x3) 

(1x3) 


If we post-multiply this matrix by the following (6x6) matrix, rank will be con- 
served. 


x 3 | u 3 I u 3 x x 3 | 0 

(3x1) | (3x1) | (3x1) , (3x3) 


U I x 3 I u 3 I U 3 X X 3 

(3x3) | (3x1) | (3x1) | (3x1) 


The product is a (4x6) matrix of rank less than 4 and, using the relations 
describing the parking orbit constraints, this matrix is: 


("AX ) 3 • x 3 


2V* 

0 

0 

(-AX) 3.G3 


(-A X ) 3 * (5 x X ) 


0 


2R 


2 


0 


0 

0 


0 

0 

0 


(-Ap) 3 *x 3 (-Ap) 3 »u 3 (-Ap) 3 *(u 3 x x 3 ) 


The vanishing of all (4x4) determinants yields the transversality conditions: 
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(a-x) 3 

• (5 3 

x i 3 ) 

= 0 


(aF ) 3 

•<=3 

X ij) 

= 0 


R 2 

<^> 3 

* U 3 

= v 2 

(AX) 3 -x 3 


The first two equations imply that the vectors ( ST)g and (Ap ) 3 are in the plane of 
u 3 and x , i.e. , the flight plane. Note also that the last matrix above has rank 3 
which confirms the requirement that our three parking orbit constraints are 
independent. 

The last two transversality conditions can be obtained by considering the 
(14x14) sub-matrix: 



We need now to join this matrix to a (6x6) one formed by considering the (5x7) 
submatrix associated with the terminal constraints and the contributions from the 
last row. Since A*** is non-singular, the (6x7) matrix, M 
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M 


(5x7) 

H 7'< 4H >4 1 X 7 T 
(lxl) . (1x3) 


has rank less than 6 . The term H 7 -AH 4 
the 49th column by adding to the 49th column the 25th column. Such column operations 
leave rank invariant. The two remaining transversality conditions are obtained from 
this (6x7) matrix. We need now to impose upon the system the desired terminal 
constraint that the vehicle be in a specified lunar orbit. The lunar orbit is defined in 
a moon-centered coordinate system and so it is convenient at this point to digress for 
a description of the system and the establishment of a few useful relations. 

We recall that the vector equation: 



(1x3) 


arises from the need to eliminate the 1 in 


*<V Y s' 


2 ,u ,v ,w ) = 
s' s s s 


I 


4 1 1 ' ^2 



0 


describes a lunar orbit which is a conic section with respect to a coordinate system 
rigidly attached to the moon and having a focus at the center of the moon. In this 
system, the position and velocity vectors of the vehicle are: 


X = 

s 

[x s' 

*s' 

’s' 

*s S 

[u s' 

V s' 

»s 1 


The coordinate system is chosen in such a manner that the conic lies in the x g y g 
plane, the major axis is co-linear with the x s axis, and the directrix associated 
with the focus at the moon crosses the positive x s axis. 

Let C be the (3x3) matrix such that the position of the_vehicle with respect 
to the moon-centered system at time t ? , the terminal time, is X g = C (x 7 ~X f ) 
where is the position of the moon with respect to the Earth-centered reference 
frame at time 1 7 . Then: 
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and 


l , - 

— C (u 7 ) + C (x 7 “X^ ) 



be a (5x3) matrix. Similarly for 



3u 3 U 

Similarly for and — - 

3u 7 3x ? 
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Now, we return to the task of extracting two transversality conditions from 
the matrix M, in which C 7 is shown explicitly. M is a (6x7) matrix which 
can be partitioned: 


dip 1 

3ijj 1 

3 \p 

3t , 

1 


(5x1) 1 

(5x3) | 

(5x3) 

f" 

— |_ 



H 7 -(AH) 4 1 

x ? T 1 

- T 

P 7 

(lxl) | 

(1x3) I 

(1x3) 


This matrix has rank less than six. Thus all ( 6 x 6 ) submatrices are singular. 
The setting of two non-trivially vanishing determinants to zero yields the remaining 
necessary conditions. But the matrix M is not in convenient form. It can, how- 
ever, be transformed into a very simple form by a sequence of multiplications by 
elementary matrices. These operations conserve rank. Before we embark on this 
procedure, the matrix M can be written in terms of the moon-centered coordinate 
system by using the relations stated above. Thus, we have: 


M = 


3 ip 
3x 


3x 


3t 


dj> 

3u 


3u I .7 

S , 3 ip 


3 1 


h 7 -(ah ) 4 


3u 


X T 
A? 


3 ip 
3x 


~r 


c + 


djp_ 

3u. 


- T 


Careful examination of the partitioned matrices will reveal which multiplications by 
elementary matrices are useful for studying the rank of M. Notice that the two 
upper right hand submatrices contain the term cbk . This suggests that postmulti- 
plication by the (7x7) non-singular matrix: g^ 


1 1 

0 

1 0 

(lxl) | 

(1x3) 

| (1x3) 

0 | 

I 

1 -c -1 c 

(3x1) 

1 

(3x3) 

| (3x3) 

0 

0 

| I 

(3x1) I 
1 

(3x3) 

j (3x3) 


will eliminate the second term in the upper right submatrix. Thus we have: 



r— 



1 


8 ij) 

3x s + 3* 

1 

3* c | 3* c 


3x s 

3t 3u 

s 

at | 
1 

8u s 1 3x s 

1 

II 


(AH) 4 

1 

1 

i 

| 

7 — “ 

- T 1 - T t T 


V 

1 

X 7 | p 7 X 7 


In a similar way, we can use the upper right submatrices to null the upper left 
submatrix. Consider the (7x7) non- singular matrix 


1 

0 

0 

(lxl) 

(1x3) | 

(1x3) 

-C _1 3{i s 

I I 

0 

at 

(3x3) 

(3x3) 

(3x1) 

i 


i 9x I 

-c" 1 s 1 
at 

1 0 1 

1 

I 

(3x1) i 

(3x3) | 

(3x3) 
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Post-multiplication yields: 


“= 1=2 = 


H* 


it c 

9u 


X T 
/ 


it c 

9x 


- T - T -1 • 
P 7 -x ? C x C 


where 


9u 


9x 


- T -1 s - T - T - 1 * -1 s 

H* = H ? -(AH) 4 -X 7 i C 1 - (p 7 X - X ? C C)C 


To eliminate the terms C, we use: 


E 3 = 


1 

0 

0 


0 

,-1 

0 c 


0 

0 

-1 


Thus: 


M*EME 1 E 2 E 3 = 


9u 


H* 


X* 


9j; 

9x 


where 


X* = (A 1 *, X 2 *, X 3 *) = X ? T C 1 


P* = (P-l*, P 2 *, P 3 *) =(p 7 T -X ? T C X C) C 
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By the way, the inverse C 
fact, will be replaced by C T 
require the specification of 


we have been using so liberally, does exist and, in 
since C is orthogonal. Further simplifications 

if . 


An arbitrary lunar orbit can be described in the following way. We choose 
our moon-centered coordinate system after the orbit is specified. We choose the 
x g , y g plane to be the plane of the orbit and so z g = 0 and r| = x| + y| . 
The constants are e, the eccentricity of the orbit, i , the semi-latus rectum and 
y^, the gravitational constant of the moon, i.e., the product of the universal 
constant of gravitation and the mass of the moon. The following equations may be 
derived directly from Kepler's laws: 



i 
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Now, 

it is useful to display 

M* completely. 






3*, 

8*, 


0 

1 

0 

0 


X 

0 






9^2 


0 

0 

1 

0 

9X 

S 


0 

0 

0 

0 

1 

0 

0 

0 





3*4 

3*4 


0 

0 

0 

0 

3x s 

3 *s 

0 

0 

0 

0 

0 

0 

0 

1 

H* 

V 

* 

<N 

V 

Pi* 

P 2 * 

p * 

P 3 


The following non-singular matrix will be used to yield an even simpler matrix. 
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=f o . Thus: 



It can be concluded here that H* = P 2 ** ~ 0 if M*E* is to have rank less than 
six. The matrix M*E* was developed from M by a sequence of multiplications 
by non-singular triangular matrices. These operations leave rank invariant. Thus 
we may direct our study of rank to the (6x7) matrix M*E* . We know its rank 
is less than six. Therefore, all (6x6) matrices are singular. Thus, H* = 0 



But replacing P 2 ** by its value from its definition, we have: 



Now, ij'p ^ 2 > and ^4 ( from the Kepler equations for the lunar orbit) are: 



The required partial derivatives are then: 
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Simplifying, we have: 



< r s +e *s> 



1*1* X S + V 


“P * Ye + P O* ( X +er ) 

+ 1 4 ~s 2 s s 

r 

s 


0 


But the equations ^ = ^4 = 0 are equivalent to: 




i 


v 

s 


y 


s 


r 


s 





- u 


s 


r + ex„ = z 
s s 

Using these relations, we have 




X * x 
1 s 


+ X 2* y s 


7 (+ Pi* 
v<t 1 


U s + p 2* V = 0 


The definitions of X * and P * are: 


x* = (x 1 * / x 2 * , x 3 *) = r ? T c " 1 

P* = (P x *, P 2 *, P 3 *) =(p 7 t -x 7 t c" 1 C) c" 1 
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and using the fact that C is an orthogonal matrix, we have: 


*“ __ rp _ fp • rp m, rp 

X _ T C x„ + (P_ T -X_ T C C) C u a = 0 


3 7 


' s 


These last two transversality conditions are then: 


h 7 - (AH 4 ) - x 7 t c t 


— T — T T * . 
Tt~ “ (p 7 _X 7 C C) 


3 x 

,T S 

' TE“ 



5 7 T C T; S 


+ 



X 7 T C T C) 


T - 
C U, 


0 
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THE SWITCHING FUNCTION 


The transversality conditions implying that the Hamiltonian is continuous at 
points 5 and 6 will now be examined closely. Thus, at point 5, we have: 

H c = H_ 

5- 5+ 

or equivalently. 


X 5- • u 5- + p 5-' 5- P 5- = 

# • 

X 5+ * U 5+ + p 5+ * x 5+ " 5+ ^5+ 


But other transversality conditions have shown that all the Lagrange multipliers are 

continuous at point 5. Since - is also continuous and , is zero, we have: 

x_ 5+ 

5 


K 5 ‘ (u 5- " u 5+ } ‘ < K 7 ) 5 P 5- = 0 


Note that — is acceleration and specific impulse is given by 


I =— . 

s p g Q e 

It was shown in Section V that: 

^ _ F 
x m 

where A = | x | 


X - 
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Then the continuity of the Hamiltonian at point 5 implies: 


Lee. 


m 


5- 


( V 5 


0 


Similarly, the continuity of the Hamiltonian at point 6 implies: 


fo 

m 


I 


6 + 




A 6 - (X 7>6 


0 . 


These conditions can be used as switching functions; points 5 and 6 can be determined 
by the vanishing of these expressions. 

If the switching function changes sign more than twice between points 4 and 7, then one 
or more thrusting periods, intermediate to points 5 and 6, may be required. 
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SUMMARY OF BOUNDARY CONDITIONS 


The conditions, which apply at the end points of the sub-arcs of an optimum 
trajectory and which cannot be trivially satisfied, are summarized below. In addition, 
the initial conditions which are unknown a priori at the endpoints of the sub-arcs are 
stated. A condition implying continuity of a quantity over a boundary point will not 
be stated and its value at the boundary point will not be considered to be free. 

Point 1 

Unspecified Quantities : ^il' ^ 31 ' ^71 ^21 is arbitrarily set equal to unity 

since all equations are homogenous in the Lagrange multipliers) . 

Point 3 


Constraints: 



(AX) 3 • (u 3 x x 3 ) = 0 
(ap) 3 * (u 3 x x 3 ) = 0 
V 2 (AX) 3 *x 3 = R 2 ( Ap ) 3 ’ u 3 


(ah) 3 = 0 


Unspecified Quantities: t 3 , X 3+ , P 3+ 
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Point 4 


Unspecified Quantity : 


Point 7 

Constraints ; ^(t^, u^, x^) =0 (i = 1,2, 3,4, 5) 
These last two transversality conditions are then: 


h ? -(ah 4 ) -a 7 t c t 


9u 

c 

W 


^ V 

- T - T T* T s 
-(p 7 T -A 7 T C T C)C T ^ 


= 0 


and 


^ g _ T 1 rn_ _ O’ _ rn m. m_ 

- -iy A ? C x g + (p ? -A ? C C)C u s = 0 


There are a total of 14 unknown quantities in the problem, precisely enough to satisfy 
the 14 conditions. Some of these conditions can be easily satisfied; for example, one 

can integrate the differential equations up until the time at which x^.x^ = R ; i.e. , t^ 

can be chosen to be that time at which the distance from the Earth's center is R. 


National Aeronautics and Space Administration 
Electronics Research Center 
Cambridge, Massachusetts, October 1966 
129-04-04-07-25 
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“The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space . The Administration 
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of information concerning its activities and the results thereof 

—National Aeronautics and Space Act op 1958 
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